Fractional S z excitation and its bound state around the 1/3 plateau of the S=l/2 

Ising-like zigzag XXZ chain. 
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We present the microscopic view for the excitations around the 1/3 plateau state of the Ising-like 
zigzag XXZ chain. We analyze the low-energy excitations around the plateau with the degenerating 
perturbation theory from the Ising limit, combined with the Bethe-form wave function. We then 
find that the domain-wall particles carrying S z = ±1/3 and its bound state of S z = ±2/3 describe 
well the low-energy excitations around the 1/3 plateau state. The formation of the bound state 
of the domain-walls clearly provides the microscopic mechanism of the cusp singularities and the 
even-odd behavior in the magnetization curve. 
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I. INTRODUCTION 

The microscopic view of the low-energy excitation in 
frustrated quantum spin systems has been on a field 
of the active researches, since the structure of the low- 
energy excitation directly reflects on the experimentally 
observable quantities such as magnetization curve. For 
instance, various interesting behaviors of the magnetiza- 
tion curve are observed experimentally for such a system 
as SrCu(B03)2 1 and the exotic nature of the low-energy 
excitation is actually reported behind the characteristic 
magnetization curved. However, in general, it is not easy 
to reveal the role of the frustration in the low-dimensional 
quantum spin systems theoretically, where the quantum 
fluctuation and frustration affects the low-energy excita- 
tion cooperatively. 

Among various frustrating spin systems, the S = 1/2 
zigzag XXZ chain is one of the most fundamental mod- 
els having the frustrating interaction without loss of the 
translational invariance^^LLS, and is actually syn- 
thesized as SrCuOziS, Cu(ampy)Brz±I, F 2 PIMNHi£ and 
(N2H5)CuCl3^. The Hamiltonian of the zigzag chain 
has a very simple structure, and has been playing a cru- 
cial role for the purpose of understanding the frustration 
effect, since it captures a variety of interesting behaviors 
induced by the frustration. In fact, the zigzag chain in 
a magnetic field has been studied actively>i^4i^iSiiLi2iiS 
Very recently, we have presented the exotic magnetic 
phase diagram of the zigzag chain including the strongly 
frustrated regionjSS we have found the magnetization 
plateau at 1/3 of the full moment accompanying the 
spontaneous breaking of the translational symmetry, the 
cusp singularities, and the interesting even-odd effect in 
the magnetization curve. 

Although the above exotic properties of the magnetiza- 
tion curve have been illustrated by the extensive numer- 
ical calculations, the microscopic mechanism for them 
around the 1/3 plateau has not been investigated sys- 
tematically. In order to address such microscopic views, 
we focus on the two key features of the zigzag chain: the 



first one is that the If J, spin structure is realized at the 
1/3 plateau state, which is extended up to the Ising limit 
of the zigzag XXZ chain22*2I. Another one is that the 
Hamiltonian of the zigzag chain interpolates between the 
single Heisenberg chain and the double Heisenberg chain 
continuously without loss of the translational invariance. 
How can we connect these two features with the exotic 
magnetization curve of the isotropic zigzag chain? A 
more thorough understanding of the problem is clearly 
desirable, which will provide an essential insight for the 
effect of the frustration in the zigzag-structured system. 

In this paper, we study the microscopic picture of the 
low-lying excitation around 1/3 plateau for the zigzag 
XXZ chain. We first classify the excitation on the 1/3 
plateau in the Ising limit and then take account of the 
quantum effect with the degenerating perturbation the- 
ory. We find that the low-energy excitation can be repre- 
sented as the combination of the domain-wall(DW) type 
excitations carrying S z = ±1/3 and ±2/3, where the 
DW particle of S z — ±2/3 can be regarded as the bound 
state of the two S z — ±1/3 particles. We further analyze 
the two-body problem of the S z — ±1/3 DWs invoking 
the Bethe type wave function^, and then clarify that the 
formation of the DW bound state becomes important as 
the double-chain nature of the system becomes dominant. 
Calculating the magnetization curve with the density- 
matrix renormalization group(DMRG) method 2 —, we fur- 
ther delineate that the DW particle picture explains well 
the low-energy excitation and the magnetization curve 
around the 1/3 plateau state. In addition to these, we 
report that the DW defect can be inserted even in the 
1/3 plateau state for a finite size system with the open 
boundary condition ( O B C ) . 

This paper is organized as follows: In the next section, 
we introduce the zigzag XXZ chain briefly. In Sec. Ill, 
we describe the 1/3 plateau state of the zigzag chain in 
the Ising limit. The basic structure of the DW excitation 
around the plateau is explained here. In Sec. IV, we take 
account of the XY term with the degenerating pertur- 
bation theory. In order to solve the eigenvalue problem 
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for the two DW problem, we invoke the Bethe type wave- 
function and clarify the bound-state formation condition. 
In Sec. V the DMRG result is presented and the relation 
with the analytical results are discussed. In Sec. VI, we 
discuss the effect of the OBC on the 1/3 plateau state. In 
Sec. VII, the conclusions are summarized. We also make 
a comment on the recently synthesized zigzag compound 
(N 2 H 5 )CuCl^. 



II. MODEL 

The model we consider here is the Ising-like XXZ chain 
with the nearest neighbor(NN) coupling Ji and the next- 
nearest neighbor(NNN) one J 2 in a magnetic field H. 
The Hamiltonian of the model is given by 

7~L = 7~to + T^zccmam (1) 

where H.q is the Hamiltonian of the system part and 
W Z eeman = —H Y^i=i &t ls the Zeeman interaction term. 
The Hamiltonian of the system part consists of the NN 
and NNN interaction terms: 7io = Wnn + T^nnn with 

3N 37V 

HnN = Jl fej.i-l-li WnNN = J2 ^ hi ; i + 2, (2) 
i=l i=l 

h itj = e(SfSf+SfS?) + SfS? (3) 

where S is the S = 1/2 spin operator, and e denotes the 
anisotropy of the XY term. In the paper, we consider the 
system of 3N sites with N =even. We also introduce the 
notation a = J2/J1 for simplicity. 

In the Ising limit (e — 0), the phase diagram at 
zero temperature was obtained by Morita and Horiguchi 
about 30 years agoi^i The 1/3 plateau appears in the re- 
gion denoted as T|| in the Fig^ The regions denoted 
as 1 | and ft II mean the Neel and double-Neel phases 
respectively, both of which have zero magnetization. An 
interesting point in the phase diagram is that the phase 
boundary of the 1/3 plateau merges the boundary of 
the Neel and double-Neel transition point at a = 1/2. 
This implies that the system is highly degenerating at 
a = 1/2. 




FIG. 1: The magnetic phase diagram of the zigzag Ising chain 
obtained in Ref. [St]. In the figure, Ji is normalized to be unity. 



For the isotropic case(e = 1), the phase diagram in- 
cluding the strongly frustrated region is recently ob- 
tained by the authors^; the 1/3 plateau survives for 
0.56 < a < 1.25 in spite of the strong quantum fluc- 
tuation. In addition, the cusp singularities show inter- 
esting behaviors around the 1/3 plateau. Moreover, the 
even-odd behavior of the magnetization curve appears 
in a large J 2 region. When analyzing these interesting 
magnetization curve, a key feature of the zigzag chain is 
that the frustration effect can be casted as the interplay 
between the single chain and double chain natures. For 
example, the two component Tomonaga-Luttinger(TL) 
liquid is realized below the low-field cusp and/or above 
the high-filed cusp, reflecting the two chain nature of the 
system, while the middle-field blanch consists of the one 
component TL liquid^ 

However the microscopic mechanism of these behaviors 
has not been illustrated yet, since the incommensurate 
nature of the system due to the frustration is difficult 
to treat theoretically. In what follows, we address the 
problem starting from the Ising limit of the zigzag XXZ 
chain. 



III. ISING LIMIT 

In this section , we investigate the nature of the exci- 
tations around the 1/3 plateau state of the zigzag XXZ 
chain in the Ising limit. Although the argument in the 
Ising limit is quite simple, we can illustrate the basic as- 
pects of the excitations around the plateau, providing a 
good starting point for taking account of the quantum 
effect with the perturbative argument of the XY term. 

A. 1/3 plateau state 

We first remark the 1/3 plateau state in the Ising limit 
briefly. The ground state phase diagram of the zigzag 
Ising chain was obtained by Morita and HoriguchiiSi As 
soon as J 2 is introduced, the 1/3 plateau appears in the 
region denoted as ||| in Fig^ Here we should recall 
that the Hamiltonian is translational invariant; at the 
1/3 plateau state, one of the following three configura- 
tion is chosen, accompanying the spontaneous symmetry 
breaking of the translation: 

iNeeb) = -ummmt- 
|Neei 2 ) = -mmmm- 
iNeei 3 ) = -mmmm-. 

For the system of 3N sites with the periodic boundary 
condition(PBC), the energy of the above Neel states is 

N 

E g = --(J 1 + J 2 ), 

where we have neglected the energy of the Zeeman term. 
The total magnetization of the state is trivially M = 
./V(or the magnetization per site is m — M/3N — 1/3). 
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type spin array S z energy 



1-2 


■■■ tittti ■■■ 


1/3 


i(Ji + J 2 ) 
|(Ji + J 2 ) 


3-1 


••• ITtUt ••• 


1/3 


3-2 


••• -LTTTT4 ••• 


2/3 


|(Ji + J 2 ) 


1-3 


■■■ 444444 ■■• 


-1/3 


|(-Ji+2J 2 ) 
|(-Ji+2J 2 ) 


2-1 


••• TTITIt ••• 


-1/3 


2-3 


•••mm--- 


-2/3 


|(Ji -2J 2 ) 



TABLE I: The low energy excitation on the 1/3 plateau state, 
which is described by the domain wall made from three type 
of the Neel ordered spin aligns. These domain-walls can be 
regarded as a quasi-particles carrying fractional S z . The num- 
bers noted in the "type" column indicate the three types of 
Neel ordered spin aligns. 



B. excitation around the plateau 

In order to analyze the low-energy excitations around 
the 1/3 plateau state, it is instructive to recall the ex- 
citation of the XXZ chain(without NNN term) at zero 
magnetic field, which is described by the DW separating 
the two-fold degenerating spin arrays of the Neel-ordered 
ground stated For the present case, the excitation on the 
plateau state is described by the DW that is made from 
the combinations of the three Neel-ordered spin aligns. 
Since the plateau state is three fold degenerating, there 
are 6(=3P2) possible types of the DW, which are sum- 
marized in TABLE I. As will be seen soon, the DW can 
be regarded as a particle carrying fractional value of S z . 
The energy of the DW particle measured from the 1/3 
plateau state is also listed in TABLE I. 

In order to see the origin of the fractional value of S z , 
we start with the wave function for M = N+l sector that 
is created by operating S + for a site in the 1/3 plateau 
state: 

u4o) = ■■■ immmmt • • • 

where "dot" is assigned for the inverted spin. The energy 
of the IV'oo) state is 

E+=E g + (J, + J 2 ). 

The \ipQ ) state is not represented as a simple DW. By 
exchanging the nearest neighboring f and J. spins, we can 
decompose the five up spin cluster into combinations of 
the fundamental DWs and thus generate the low-energy 
excited state for the M = N + 1 sector systematically: 



l^di) = ■■■imttmtmmt--- 

l^o + 2 > = ■•• 44444444444444444 •• • 
|V>i + 2 ) = • ■ • 44444444444444444 • • • 

3/ sites 3m sites 

\i>L) = ■■■ T4TTT4 • ■ • 444 444444^444444 • • • (4) 



where we assign underlines to the | i T or T I T clusters 
that are useful signs to find the DWs. For all configu- 
rations above, we can easily confirm that the energy is 
the same as E + = E g + (Ji + J 2 ). These degenerating 
states are connected by of XY term, which is essential 
for the degenerating perturbation treatment in the next 
subsection. 

By comparing the spin configurations in eq. J3J with 
those in TABLE I, we can identify the 444 cluster as the 
DW carrying S z — 1/3; Since all of the 444 clusters in 
l^lm) s t a te are equivalent, the three DWs are sharing 
the magnetization S z — 1 equally. Thus each DW is 
represented as a particle of S z = 1/3 and AE = (Ji + 
J 2 )/3m^ In addition to this, the configuration of 4444 
can be regarded as a bound state of the two S z = 1/3 
DWs, which has S z = 2/3 and AE = 2(Ji + J 2 )/3. As is 
described in the next section, this bound state becomes 
important for J 2 3> J\. 

Also for M = N — 1 sector, we can make almost the 
parallel argument to the S z — 1/3 DW case. Inverting 
an up spin in the 1/3 plateau state, we have a state: 

|^oo> = ••■ 44444444444444444444 

whose energy is E~ = E g . Although this value of the en- 
ergy is the same as that of the 1/3 plateau state, the total 
energy including the Zeeman terms for the 1/3 plateau 
state is lower than that for I^qq) in the certain range of 
the magnetic field. 

We also generate the low-energy states for the M = 
N — 1 sector, by exchanging the nearest neighboring 4 
and 4 spins. However, we note that there is a bit different 
point from the M = N + 1 case. Since the state \ip^ ) 
contains two kind of the spin clusters, 44 an d 444 > wc 
introduce the slightly modified labeling for the position 
of the 44 cluster; we assign "bar" to the label of the 44 
configuration. 

I^oi) - •■■ 44444444444444444444 ■• ■ 

l^oi) = ••• 44444444444444444444 •• • 

l^u) = ••• 44444444444444144444 ■• ■ 

Ife) = •••44444444444444444444 ••• 

3/ sites 3m sites 

HrJ = •■■44444^^444444^444444 ••• 

where underlines are assigned for the 444 an d 44 clusters. 
Clearly the 444 cluster can be regarded as a DW-particle 
carrying S z = —1/3, corresponding to the type 1-3 or 2-1 
DWs in TABLE 1. Moreover, the energy for the states 
consisting of the three 444 DWs is easily calculated as 
E = E g - Ji + 2J 2 . Thus the energy of the S z = -1/3 
DW particle reads AE = (- Ji + 2J 2 )/3. 

We further consider the 44 cluster, which is depicted 
as the 2-3 type DW in TABLE 1. Since the 44 cluster is 
decomposed into the two 444 clusters, the 44 cluster can 
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be regarded as the bound state of the two S z = —1/3 
DWs. Thus we find that the J, J, cluster is the DW bound 
state carrying S z = —2/3. In addition, we note that 
the total energy for the state including the S z = —2/3 
DW(for example IV'm.)) ^ E = E g . Therefore the energy 
of the S z — —2/3 DW bound state is readily identified 
to be AE = (Ji - 2J 2 )/3. 

An essential point on the 1 /3 plateau state is that the 
low-energy excitations of the M = N ± 1 sector are de- 
scribed by the combinations of the DWs listed in TABLE 
1. We next take account of the quantum effect originat- 
ing from the XY term on the basis of the DW picture. 



As was seen above, if the NNN term is operated on 
the single DW state, it generates only higher energy con- 
figurations. Then a natural question arises: what is the 
role of the NNN term? In order to see it, we examine 
the effect of the NNN term on the DW bound state of 
S z = 2/3. We define the state for the DW bound state 
at the position z by 



\z) 



mmmm 



(8) 



Then the matrix element of 7Ynnn for the DW bound 
state is 



IV. PERTURBATION BY THE XY-TERM 

In this section, we consider the the quantum effect on 
the DW excitations illustrated in the Ising limit, by using 
the degenerating perturbation theory with respect to the 
XY term. We assume e <C 1, for which the perturbation 
treatment is justified. Here we should note that the cor- 
rection to the energy of the 1 /3 plateau state starts with 
2nd order of s: 



E g /N 



1 



( Ji + Ji 



~2 



J l 
J 2 



Ji 



(5) 



Thus we can set the origin of the energy as E g — 
— + ^2) in the following arguments of the first order 
degenerating perturbation for the DW excitations. 



A. basic views for the DW 

We analyze the one-body problem of S z = 1/3 DW 
in the infinite length chain up to the first order of e. 
Although the higher order terms may be required for a 
quantitative analysis of the problem, we can capture the 
intrinsic property of the DW excitation within the first 
order theory. 

We label the position of the S z = 1/3 DW as x which 
is defined as the site of the center of three up spins ^& 
Then the matrix element of the 7io is obtained as 

n \x) = ^±^\x) + £ -^(\x-3) + \x + 3)) 



Hnnn^) = |(Ji + J 2 )\ z ) + ^(\z-3) + \z + 3)) 

+higher energy terms (9) 



where the higher energy configurations are generated by 
the NNN term operated on the spins away from the DW 
part. Thus we see that the DW bound state can move 
by using the NNN term without the extra energy cost. 
This is a crucial point on the NNN term. When J\ ^> J2, 
the low-energy excitation on the plateau is basically de- 
scribed by the single DW excitation. On the other hand, 
when Ji <C J2, we can expect that the role of the DW 
bound state becomes essential in the low-energy excita- 
tion. The switching mechanism between the single DW 
excitation and the DW bound state is clearly associated 
with the crossover between the single chain nature and 
double chain nature of the zigzag chain. In the next sub- 
section, we analyze the two DWs and its bound state 
problem systematically. 



B. two body problem and bound state 

In order to understand the formation mechanism of the 
DW bound state, we consider the two body problem of 
the DWs in the infinite chaim 22 Introducing the notation 
for the free two DW state, 



+ higher energy terms 



(6) 



\x,y) 



mm? •••mm? 



where the higher energy terms include the configurations 
that yield the energy rise of order J\ or J 2 . Here an 
important point is that the NNN term always gener- 
ates higher energy terms. Thus the NNN term can not 
contribute to the low-energy dynamics of the single DW 
within the first order of e. By neglecting the higher en- 
ergy terms, we can easily obtain the dispersion curve of 
the single DW: 

uj(k) = l(J 1 + J 2 )+eJ 1 cos{3k). (7) 
o 



we write the wave function as 



W) = ^2fx, y \x,y) + ^g z \z) 



x<y 



(10) 



where the summentions about x,y,z are taken over the 
all possible positions of the DWs. Then we obtain the 
real-space eigenvalue equation for the two DW scattering 
problem: 
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2 eJ% 

E^fx,x+4 = g (Jl + Jl)fx,x+& H ^- [/a;-3 : x+4 + fx,x+7 + .9k + 5^+3] j 

E {2) g x = ~(Ji + <h)9x + [9x+3 + 9x-3] + ^ [fx-3,x+i + fx,x+4 ■ 

In order to deal with eqs. and l(12|). we assume the Bethe type wave function for the two DW problem: 

f x>y = A(k, k')e lkx+lk ' v + B{k, k')e lky+lk ' x . 
k,k' 



We also write the wave function for the bound state as 



9x = U ^ k ') 



^i(k+k )x 



The energy eigenvalue in eqs. 1(11(1 and (|12|l is given by 

£ (2) (fc,fc') =uj(k)+uj{k') 1 



(11) 
(12) 

(13) 

(14) 

(15) 



which is determined from the free part of the two DW problem. 

We define the S'-matrix of the two DW as S(k, k') = B/A. Solving the scattering problem 111(1 and 112(1 . we obtain 



S(k,k') = -e l ( fc '- fc )- 



g — i3k _|_ ^i3k' 

k) 


2(cos(3fc) - 


hcos(3/c'))-a(e i3(fe+fc ' ) 


_|_ g-i3(fe+fc')) 


£j 1 3 j g i- 3 A, 


2(cos(3fc) - 


h cos(3fc')) - a(e i3 ( fc + fc ') 


|. e -i3(fc+fc')y 



(16) 



r 



In this expression of the S'-matrix, the extra phase factor 
e i(k -fc) means if a overtakes another one, the 
spin array around the DW is shifted by one site. This is 
a peculiar factor for the DW problem, unlike to the usual 
spin wave one. 

We now consider the formation of the bound state in 
the two DW problem; we assume 



k = u + iv and k! = u — 



(17) 



which corresponds to the 2-string solution in the Bethe 
ansatz terminology. Substituting eq. 1(17(1 into eq. 1(15(1 , we 
obtain the dispersion curve of the bound state particle 



E 



bound 



(«) = o ( J i + J 2) + EJi(e~ iv + e 6v ) cos3u. (18) 



From the normalizability condition of the wave func- 
tion, u and v must satisfy a nontrivial condition: 



e 6v cos(3u) + 2acos^(3u) - a = 0, 



(19) 



which is equivalent to the pole of the S- matrix 1)16(1 . The 
formal solution of eq. 1(19(1 is given by 



„3u 



cos 6u 
* 

COS 6U 



(20) 



Since the physical solution must satisfy e 3v > 1 and 
— 1 < cos3w < 1, the bound state dispersion curve 
emerges in the restricted range of u. We numerically eval- 
uate ea. l(20|) in the permitted range of u with e = 0.1, and 



illustrate the obtained dispersion curve l(18|) in FIG [21 In 
the figure, we can see that the bound state band appears 
below the dispersion curve of two free DWs that is ob- 
tained as eq. ((15(1 with fc = fc'. 

On the bases of figure |3 we discuss features of 
E hound (u) in detail. The low-energy solution of ea. l|2l7|) 
exists in the range 

- < u < u or u <u<~, (21) 

6 3 2 K ' 

where u* = | arccos[| (a -1 — y 'a~ 2 + 8)]. In addition 
to this, we note that, when a > 1, ea. 1(19(1 has a solution 
in the high energy region around u ~ and 2n/3. 

In FIGI2 we can see that, as a increases, the bound 
state dispersion curve comes down to the lower-energy 
region. Since the minimum energy of the free two-particle 
dispersion JTSJl is given by E^ n = 2(Ji + J 2 )/3 - 2eJi 
at fc = fc' = 7r/3, the bound state excitation becomes 
relevant, if the minimum of the dispersion ((18(1 is lower 

than E£( n . Substituting eq.® into eq.©, we obtain 
the range of a where the bound state is relevant in the 
low-energy region: 



a > 2 



al). 



(22) 



This value of a+ is particularly important for the even- 
odd behavior of the magnetization curve, which is dis- 
cussed in Sec. V. 

In addition, we note that the minimum energy of 
the DW bound state is achieved at u ~ 7r/6 or w/2, 
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J 2 =1.0 




1/3 71 u 2/3 jc 



FIG. 2: The DW bound state dispersion curve for e — 0.1 
and Ji = 1. The curves in the low-energy region are shown 
for < u < 2n/3. The solid lines indicate the curves for 
Ji =1.0, 0.8, 0.5(= l/a+) and 0.4. The dashed lines are the 
corresponding lower bounds of the free two DWs. 

where cos3u ~ 0. Thus we can easily see E hound (u) ~ 
2(Ji + J2)/3 + J 2 cos 6m near the bottom of the bound- 



state dispersion curve. This implies that the low-energy 
excitation for J 2 ~> J\ is described by the hopping of the 
DW bound state originating from the NNN term, which 
is consistent with the basic view of the DW bound state 
discussed at eq.©. 



C. two body problem for the S z = -1/3 DW 

We next consider the two-body problem for the S z — 
— 1/3 DWs. Within the first order perturbation, we can 
solve the two DW problem with the same manner as the 
previous S z = 1/3 case. 

The single DW dispersion curve is obtained as 

uj(k) = 1 ^ 2 + £ J\ cos 3fc (23) 

The S'-matrix of the two-DW problem is also calculated 
as 



S(k,k') 





2A + 2(cos(3fc) - 


hcos(3fc'))-a(e l3( ' £+fe ' ) - 


_ e -i3(k+k')\ 


^j23& j g 2.3/^ 


2A + 2(cos(3fc) - 


hcos(3fc')) - a(e l3 ( k + k ') - 


i_ e -»3(fc+fc')y 



(24) 



where 



-l + 2a 



(25) 



This A terms originates form the "binding energy" in the 
zcroth order(Ising limit). 

Assuming the 2-string type solution l(T7|) . we obtain 
the equation for the DW bound state: 



where we can actually see that the behavior of the curve 
changes drastically. 

When a < |5f > the solution of ea. H26l) appears only in 
the high energy region(it ~ 0). As a increases, the bound 
state excitation is able to appear in the low-energy region; 
For 



a > 



1 



(29) 



e Ju cos3u + 2a;cos 2 3w-a- A = 0, (26) 



and its formal solution: 



a cos 6u — A 
cos3u 



(27) 



The dispersion curve of the DW bound state is 



E hound {u) = -(-J 1 + 2J 2 )+eJ 1 (e- 3l, -re 3, ')cos3u. (28) 

Although this expression of the dispersion looks very sim- 
ilar to the S z — 2/3 case, the resulting behavior of the 
curve is quite different; Since the A contains the factor 
l/e, the dominant properties of (|2*7Jl is well approximated 
to be e 3v ~ —A/ cos 3m as far as |A| > e. However, for 
A < e, namely, in the vicinity of a — 1/2, the curve 
depends on a sensitively. We thus show the bound state 
dispersion of e = 0.1 around a = 1/2 in FIGE3 



the bound state dispersion curve appears at u = 7r/6 
and 7r/2. As a is increased beyond this value, the two 
branches of the curve are extended from the band edges 
toward u = 7r/3. At the same time, the energy at u = 7r/6 
and 7r/2 decreases rapidly. We can see that E hound (u) has 
the local minimum at the band edges for 



a > 



1 + y/1 + 2e(2 + e) . 



2(2 + e) 



(= a-), 



d 2 £ bound ( M ) 1 
du 2 



\u= 7r/6 or 7r/2 



(30) 



= 0. 



where a_ is determined by 

Since the minimum value of the two free DW dispersion 

(2) 

is given by -E„ lin = 2(— J\ +2 J2)/3— 2eJ\, we can further 
see that the bottom of the bound-state dispersion curve 



(2) 

becomes lower than E^J n , when 



a > 



l + 2e 



(= «-)■ 



(31) 
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J 1= 1.0 




2/3 Jt 



FIG. 3: The DW bound state dispersion curve of e = 0.1 
around a — 1/2. The curves in the low-energy region are 
shown for < u < 2tv/3. The solid lines indicate the curves 
for J 2 =0.5, 0.53, 0.55, 0.571(~ at), 0.58 and 0.6. The dashed 
line is the lower bound of the two free DW dispersion. In this 
figure, the origin of the energy axis is shifted with £; boulld — 
E h 



+ 2( Ji - 2J 2 )/3. 



This value of a*_ is also important for the analysis of the 
even-odd behavior of the magnetization curve. 

As a is increased further, after the branches extended 
from u = 7r/6 and n/2 connect with each other at 
a = |^|, E bound (u) approaches ±(Ji - 2J 2 ) + eJ 2 cos6u 
rapidly. This behavior of E hound (u) is also consistent 
with the DW bound state moving with use of the NNN 
term when J-x 3> J\, 



V. MAGNETIZATION CURVE 

In order to verify the analytical results for the DW ex- 
citations around the 1/3 plateau, we calculate the mag- 
netization curve for the Ising-like zigzag XXZ chain with 
the DMRG method. We discuss the relevance of the 
fractional S z excitations to the obtained magnetization 
curves. 



A. numerical results 

We calculate the magnetization curve for the zigzag 
XXZ chain of 192 sites(7V = 64) for e = 0.1, using 
the DMRG method. The number of the retained bases 
used in the DMRG computations is typically 64, with 
which the calculated magnetization curves converges suf- 
ficiently. 

Let us first survey the features of the curves. In FIG0J 
we show the magnetization curves for J± > J2, where we 
fix J 1 — 1 and vary J 2 . As J 2 is increased from J 2 = 0, 
the magnetization plateau emerges at 1 /3 of the full mo- 
ment(FIGlU-(a)), and the width of the 1/3 plateau ex- 
tends rapidly. At the same time, the gap at zero magnetic 
field becomes small. This is consistent with the transition 
form the Neel(||) phase to the double-Neel(TTII) one in 
the Ising limit, which is located at a = 1/2 and H = 0. A 



precise study of the critical point of the Neel and double- 
Neel transition for the zigzag XXZ model can be seen 
in Ref, 2 — . The estimated value of the critical point for 
e = 0.1 is J2 ~ 0.45, which is in good agreement with 
the magnetization curve in FIG0J(b). In addition to this, 
we can see that the low-field branch(= the magnetization 
curve below the 1/3 plateau) becomes very narrow. As 
was mentioned in Sec. II, the zigzag Ising chain is highly 
degenerating at a = 1/2 below the 1 /3 plateau, where the 
phase boundary of the 1 /3 plateau also merges with the 
Neel and double-Neel transition point (see FIG^I. For 
the XXZ chain of e 7^ 0, the quantum fluctuation com- 
ing from the XY term lifts this degeneracy, and thus the 
magnetization curve appears in the narrow region be- 
tween the zero magnetization and the 1/3 plateau. The 
detailed relation of the DW excitation and the magneti- 
zation curve is discussed in the next subsection. 

Turning to the high-field branch of the magnetization 
curve, we can see that the cusp singularity appears in 
the high-field branch for J 2 > 1/4, which can be ex- 
plained well by the shape change of the spin-wave dis- 
persion curve from the saturated stated 

In FIGQ}(c), we can find that the even-odd oscillation 
appears in the low-field branch. Although the true mag- 
netization curve of the finite size system should be deter- 
mined from the stable M =cvcn states, we here show the 
anomalous M =odd steps as well, because it is the evi- 
dence of the DW bound state having S z = ±2/3. We dis- 
cuss the detailed mechanism of the even-odd oscillation 
based on the DW picture in the next subsection. Here we 
note that the even-odd behavior emerges for J 2 > 0.58. 
As J2 is increased further, the amplitude of the even- 
odd oscillation in the low-field branch is extended grad- 
ually, and at the same time the cusp singularity in the 
high-field branch shifts to the low magnetization side (see 
FIGH(d)). 

We further show the magnetization curves for J\ < J 2 
in FIG El In the figures we fix J 2 = 1 and vary J\. As 
Ji is decreased from unity, the width of the 1/3 plateau 
shrinks. Simultaneously the region of the even-odd oscil- 
lation in the low-field branch extends down to the zero 
magnetization, and the amplitude of the oscillation be- 
comes significant (FIG 0- (a)). On the other hand, for the 
high field branch, the position of the high-field cusp ap- 
proaches the 1/3 plateau; Indeed for the magnetization 
curve of J\ ~ 0.6 in FIG|Sl-(b) we can see that the high- 
field cusp merges into the 1/3 plateau. For J\ < 0.6, 
the cusp singularity does not appear, but, instead, the 
even-odd oscillation emerges in the high-field branch. In 
FIG [5J (c), we can confirm the oscillating behavior of the 
magnetization curve above the 1/3 plateau clearly. In 
Ji — > limit, the 1/3 plateau disappears and the magne- 
tization curve finally becomes that of the two decoupled 
XXZ spin chain. Here it should be remarked that, in the 
high-field region, the phase diagram of the Ising limit 
contains no special point having a peculiar degeneracy, 
unlike to the previous Neel and double-Neel transition 
point in the low-field case. This suggests that the even- 
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FIG. 4: The magnetization curves of the Ising-like zigzag XXZ 
chain of the 192 sites for e = 0.1. (a)J 2 = 0.2, (b)0.5, (c)0.6 
and (d)1.0. The inset in figure (c) is the magnification around 
H ~ 0.2, where the even-odd oscillation appears in the low- 
field branch. 

odd oscillation in the high-field branch is induced by the 
purely quantum effect. 

B. DW excitation and the magnetization curve 

We now discuss the relevance of the DW excitations 
to the magnetization curves. The magnetization curve 




1 H 2 



FIG. 5: The magnetization curves of the Ising-like zigzag XXZ 
chain of the 192 sites for e = 0.1. (a)Ji = 0.8, (b)0.6, and 
(c)0.4. The inset in figure (c) is the magnification around 
H ~ 1.3, where the even-odd oscillation appears in the high- 
field branch. 



of the ID quantum spin systems is generally described 
by the hard-core bosonic particle picture^ 28,29 a particle 
having magnetization fills its dispersion curve up to the 
"chemical potential" corresponding to the external mag- 
netic field. Then the magnetization curve is interpreted 
as the chemical potential v.s. the particle-number curve, 
where the shape of the dispersion curve are essentially 
important to figure out the feature of the magnetization 
curve. As was seen in Sec. Ill, the elementally excita- 
tions around the 1/3 plateau are described by the DW 
particles carrying S z — ±1/3 and their bound states hav- 
ing S z = ±2/3, although the total magnetization of the 
system always takes an integer value particularly for a 
finite size system. This implies that the magnetic excita- 
tion of S z = ±1 around the plateau state is described by 
the combinations of the fractional S z DW excitations. 

Let us first discuss the relation between the S z = +1/3 
excitation and the characteristic properties of the high- 
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free DW dispersion 




FIG. 6: Schematic diagram (a) for the cusp singularity and 
(b) for the even-odd oscillation of the magnetization curve. 
The open circle denotes the DW excitation and the solid circle 
denotes the DW bound state. The broken lines indicate the 
"chemical potential" corresponding to the magnetic field 

field branch. If Ji > J 2 , the effect of NNN term is 
not so big. Thus the low-energy excitation around the 
plateau state is basically described by the free S z = 1/3 
DW excitation. The magnetization curve near the 1/3 
plateau is reflecting the shape of the u(k) for the S z — 
1/3 DW; the magnetization curve raises from the 1/3 
plateau with the square-root behavior associated with 
the curvature around the bottom of cj(k)£& This is the 
case for J 2 < 1/4 in the high field branch(FIG0Ka)) 

As J2 is increased and the NNN effect becomes more 
significant, the dispersion curve of the DW bound state 
comes down to the low-energy region. Then there are 
two possible situations: the bottom of the DW-bound- 
state dispersion curve is slightly higher than that of the 
single DW dispersion curve(FIGEK a ))j and the opposite 
(FIGEKb)) 

For the former case, the magnetization increases along 
the single DW dispersion curve, as long as the chem- 
ical potential is below the bottom of the bound state 
dispersion curve. However, when the chemical potential 
touches the bottom of the bound-state dispersion curve, 
the magnetization curve captures the band edge singu- 
larity to have the cusp (FIG 0- (a)). After the chemical 
potential exceeds the bottom of the bound-state disper- 
sion curve, the magnetization can increases rapidly by 
using the bound-state dispersion curve. In the context 
of the crossing points of the chemical potential and the 
dispersion curves, the above obserbation implies that the 
two component TL liquid is realized above the cusp sin- 
gularity, while the one component TL liquid is realized 
between the high-field cusp and the 1/3 plateau. These 
behaviors of the each branch of the magnetization curve 
is consistent with the known results based on the shape 
change of the spin-wave dispersion curve from the sat- 
uration limiU^. The situation illustrated in FIG|6j-(a) 
corresponds to the case for a < a*,. 

For the latter case, the magnetization curve increases 
by using the bound state dispersion curve. We here recall 
that the excitation of S z = 1 is always represented as a 
combination of three DW particles. If M =odd, there 
are 3(M — N) number of S z = 1/3 DW particles in the 
system; 3 -. Then 3(M — N) — 1 numbers of the DWs can 
conform the bound state of S z = 2/3, but the remaining 



one DW can not find its partner. Thus the remaining one 
DW have to sit on the free DW dispersion curve, as in 
the diagram of FIGEIO 3 )- Therefore the M =odd state 
has a slightly higher energy due to the gap between the 
bottom of the single DW band and the chemical potential 
lying in the bound state band. On the other hand, for the 
case of M =even, all of the DWs can find their partners 
and conform the bound states successfully. Clearly this is 
the origin of the even-odd behavior of the magnetization 
curve. 

According to the analysis in Sec. IV, we can see that 
the switching of the above two cases occurs at a — a* + = 
2. Although the theoretical prediction oT l = 1/2 is con- 
sistent with the DMRG result a -1 ~ 0.6, they shows a 
slight deviation. The reason for this is that the S z = 1/3 
DWs have no "binding energy" in the zeroth order(Ising 
limit), corresponding to A = in eq. (TBJ. Thus the 
leading term for the binding energy is of first order of e, 
implying that the formation of the DW bound state is 
responsible for the purely quantum effect. For the pre- 
cise determination of a* +1 the second order calculation 
with respect to e is required, where the application of 
the Bethe form wave function is rather difficult. 

We turn to the analysis of the low-field branch, which 
can be explained by almost the same line of the argu- 
ment. However, in contrast to the high-field branch, the 
S z = —1/3 DW excitation has the zeroth order binding 
energy, namely A term, which yields the sensitive be- 
haviors of the dispersion curve of the DW bound state 
near a = 1/2. Thus we first consider the case where 
a is sufficiently away from 1/2, where the shape of the 
bound state dispersion is rather simple. When a < 5x51 
the bound state DW does not appear in the low energy 
region. Thus the magnetization curve can be described 
by the free S z = —1/3 dispersion curve(FIGQJ(a) and 
(b)). On the other hand, if J2 > 0.62, the shape of 
the bound state dispersion curve is approximated well 
by ~ jj(—Ji + 2J2) + eJ2COs6fc, whose bottom is 
lower than that of the single DW dispersion. Thus the 
low-field branch of the magnetization curve exhibits the 
even-odd oscillation for J 2 > 0.62(FIGi[d) and FIG©, 
due to the same mechanism in FIGEKb). 

We next proceed to the analysis for a ~ 1/2. As is 
seen in FIG El a characteristic feature of the S z = —2/3 
DW bound state is that the shape of the dispersion curve 
changes drastically in narrow range of a. Accordingly, we 
can expect that the shape of the low-field branch of the 
magnetization curve also changes in the corresponding 
region of a. In order to verify this expectation, we calcu- 
late the the low-field branch of the magnetization curves 
for 0.5 < a < 0.62 intensively, which are shown in FIG0 
When a < a_ = 0.521, the small dispersion curve of the 
bound state appears accompanying a local maximum at 
the band edges(w = 7r/6 or u = 7r/2), which is located 
in a sufficiently higher energy region than the bottom of 
the free DW dispersion curve. Thus the magnetization 
below the 1 /3 plateau decreases using the free DW band. 
For &- < a < a*_ = 0.571, the bound state dispersion 
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FIG. 7: The low-field branch of the zigzag XXZ chain near 
a = 1/2: (a) a = 0.53, 0.55 and 0.57, and (b) a = 0.59 and 
0.62. 



comes down to the low energy region and has the local 
minimum at the band edges. Nevertheless the energy at 
this local minimum is still higher than the bottom of the 
free DW dispersion. Thus, in this region of a, the sit- 
uation is the same as the case of FIG|Hl-(a), so that the 
magnetization curve has the cusp singularity (FIG [3(a)). 

For a > a*_, the lowest energy of the bound state dis- 
persion becomes lower than that of the single DW one, 
and thus the magnetization curve exhibits the even-odd 
oscillation, which is the case of FIGEKb) and FIG[3(b). 
However, we should notice that the bound state disper- 
sion curves has another feature especially for 0.579 < 
a < 0.605;— the bound state dispersion curve has the lo- 
cal minimum at u = 7r/3, implying that the curve shows 
the double-well shape. Thus the magnetization curve 
may have the cusp singularity and even-odd behavior at 
the same time; the DMRG result for a = 0.59 and 0.6 
actually shows the both of the cusp and even-odd os- 
cillation(FIG|3(b), see also the inset of FIGH(c)). For 
a > 0.605, the bound state dispersion has the single min- 
imum at u = 7r/6 or 7r/2, and then the magnetization 
curve exhibits the even-odd oscillation without the cusp. 

As seen above, the DW excitation picture of fractional 
S z can explain well the magnetization curve calculated 
by the DMRG. In the remaining part of this section, we 
shall discuss to what extent the present theory can re- 
tain its validity. Since the dispersion curves in Sec. IV are 
obtained within the first order perturbation, the higher 
order contribution with respect to e may affect the shape 
of the dispersion curves. However, what we want to em- 
phasize here is that the present DMRG results actually 
justifies the DW picture as far as e = 0.1. Moreover, it 
should be noted that, even for the isotropic zigzag chain, 



the topology of the magnetic phase diagram is almost 
the same as the present Ising-like XXZ model, except for 
the dimer-gapped phase at the zero magnetic field of the 
isotropic case. The fact of the 1/3 plateau existing for 
the isotropic zigzag chain suggests that the DW picture 
is basically maintained against the quantum fluctuation 
due to the XY term. Of course, a further precise analy- 
sis is required for the quantitative understanding of the 
problem especially around a ~ 1/2, which may be an 
interesting future subject. 

Another important factor which should be added in 
the present theory is the interaction effect between the 
DWs. Although the DMRG results also demonstrate 
that the DW picture basically holds up to the certain 
range of the magnetization curve, a careful consideration 
of the interaction between DWs is required for a more 
quantitative analysis for the magnetization curve away 
from the 1/3 plateau. For example, the DWs and their 
bound state picture for the even-odd oscillation suggests 
that the magnetization curve may have some singularity 
at the upper edges of the even-odd oscillation branch, 
where the chemical potential touches the bottom of the 
free DW dispersion curve. However we can not confirm 
such singularity within the DMRG result for the 192 sites 
chain. A possible reason for this is the interaction effect, 
since the number of the S z = 1/3 DWs increases as the 
magnetization is increased away from the 1/3 plateau. 



VI. OBC AND DEFECT IN THE 1/3 PLATEAU 
STATE 

A. Ising limit 

In this section we give a remark about the plateau state 
for the OBC system. Since the DMRG gives the results 
for the OBC, the plateau state for the OBC has a special 
importance for the alanysis of the DMRG results. 

For the case of the PBC, it is not permitted to insert 
single DW into the plateau state without changing the 
total magnetization. In other words, the low-energy ex- 
citation of the M — N sector is always represented as a 
combination of two DWs such as the pair-creation of the 
S z — ±1/3 particles. Thus the excitation that does not 
change the total S z of the system always accompanies 
the energy rise of order J\ or Ji . 

For the case of the OBC, however, we can insert single 
DW into the plateau state with keeping the magnetiza- 
tion of the system. Thus we should check the energy of 
not only the "uniform Neel state" but also the state in- 
serted the single DW. The three uniform Neel states for 
the OBC are: 

|nmi) = mm ••■mm 
|Neei 2 ) = mm •••mm 
|Niei 3 > = mm ■•■mm. 
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and the corresponding energies are easily obtained as 



Ei — E q i + — 



= — E„ 



Since the DWs of S z = —1/3 or —2/3 have "negative 
energy" depending on the coupling constants, we should 
evaluate the energy for the following three states: 



is) = nt • 


■nun- 


•ITT 


21) = ITT • 


■mm- 


■TIT 


23) = TTI • 


■TTUTT- 


•ITT, 



which may have lower energies than those of the uniform 
Neel states <|32l The energies for these states are easily 
calculated as 



E, 



E- 



23 



E 2 i — E, 



E„ 



Ji 
9 4 

Hi- -Ii 

4 2 ' 



(33) 
(34) 



Thus we can conclude that the 1/3 plateau state for OBC 

is |Neeli) for J 2 < J\ and ^23) f° r Ji > J\- Although 
this is the case of the Ising limit, we will see that the 1/3 
plateau state that is inserted the DW survives against 
the quantum fluctuation. 



(a) 4-1.0, J 2 =0.8, e-0.1, M=32 




FIG. 8: The distribution of the local spin moment m(i) at 
the 1/3 plateau state: (a) Ji = 1.0 and J2 = 0.8, and (b) 
Ji = 0.8 and J 2 = 1.0. 



VII. SUMMARY AND DISCUSSION 



B. DMRG result 

According to the plateau quantization condition, 
the period of the 1/3 plateau state must be q = 
3 x integer with the spontaneous symmetry breaking of 
the translation^ The three kind of the Neel states clearly 
satisfy this condition. However for the OBC system, we 
have discussed that the S z — —2/3 DW can be inserted 
in the 1/3 plateau state in the Ising limit. In order to 
investigate the effect of the XY term on the defect, we 
calculate the distribution of the local moment along the 
chain with the DMRG method, which is shown in FIG[H] 
We can clearly see that the defect is inserted at the cen- 
ter of the chain(FIG|HKb)). A precise computation yields 
that the defect appears for a > 0.923, which is consis- 
tent with the result for the Ising limit. Moreover, we note 
that the DW inserted plateau state survives even for the 
isotropic zigzag chain, for which the appearing point of 
the DW-inserted plateau state is a ~ 0.85. 

Here it should be remarked that the translational sym- 
metry of period 3 for the 1/3 plateau is broken by the 
inserted DW defect. However, this does not mean the 
break down of the plateau quantization condition that is 
valid for the infinite length chain. The present DW defect 
in the plateau is the single impurity problem in a finite 
size system, which can be neglected in the bulk limit. 
Therefore, the DW defect in the plateau state does not 
conflict with the quantization theorem which is relevant 
for the infinite length chain. 



In this paper, we have studied the excitation around 
the 1/3 plateau state of the 5=1/2 Ising-like zigzag 
XXZ chain. We have explicitly constructed the low- 
energy excitations as the DWs separating the three-fold 
degenerating Neel-ordered spin arrays of the 1/3 plateau. 
The important point is that the DWs are regarded as the 
quasi-particles carrying the fractional value of S z : the 
S z = ±1/3 DW and its bound state having S z = ±2/3. 
Further, we have taken account of the quantum effect 
of the XY term with the degenerating perturbation the- 
ory. Then wc have found that the low-energy dynam- 
ics of the DWs is significantly affected by the interplay 
between the single and double chain natures of the sys- 
tem originating from the zigzag structure of the model; 
The S z = ±1/3 DW moves by using the NN interaction 
term, while S z = ±2/3 DW bound state uses the NNN 
one. Analyzing the formation mechanism of the DW 
bound state precisely, we have shown that the crossover 
of the dispersion curves in the low-energy region occurs 
at a — a* + for the S z = 1/3 DW, and at a = a*_ for the 
S z — — 1/3 DW. Moreover, this crossover of the disper- 
sion curves certainly induces the exotic behavior of the 
magnetization curve: The cusp singularity and the even- 
odd oscillation of the magnetization curve are switched 
with each other at a = a* ± , depending on the relative 
position of the dispersion curves of the free S z = ±1/3 
DW and the S z = ±2/3 DW bound state. In particular, 
the fact that the usual S z — ±1 excitation is represented 
as a combination of the S z — ±1/3 DW particles is es- 
sential for the even-odd oscillation of the curve. In fact 
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the DW excitation combined with the hard-core bosonic 
particle picture has explained well the DMRG calculated 
magnetization curves. Although the present argument 
is based on the degenerating perturbation theory of the 
first order with respect to e, the intrinsic properties of the 
DW excitations are vividly illustrated within the present 
theory, which is clearly supported by the DMRG results 
for the magnetization curve. 

Here it is worthwhile to remark about the connection 
of the present DW picture to other view points for the 
magnetization curve. As was seen above, the fractional 
S z DW is valid for the magnetization curve near the 1/3 
plateau. On the other hand, we know that a good picture 
to describe the magnetization curve in high-field region is 
the spin wave excitation from the saturation field, which 
also explains the cusp singularity successfully^. Then 
both of the DW picture around the 1 /3 plateau and the 
spin wave picture from the saturated state seem to cap- 
ture the qualitative feature of the magnetization curve 
in the crossover region of these two pictures. A similar 
situation is also considered for the low-field branch. For 
e = 0.1 the ground state at the zero magnetic field is the 
Neel or double-Neel order. More interestingly, the dimer- 
singlet order is developed particularly around a = 1/2 as 
the system approaches isotropic point (e — > 1). For the 
isotropic system, the excitation at the zero field is repre- 
sented as the DW separating the two-fold degenerating 
dimer-single states, for which the bound state of such 
DW is also reportedM^. However the relationship be- 
tween the DW picture at the 1/3 plateau and the DW 
based on other pictures has not been clarified yet. The 
unified view for these complemental pictures may be a 
key problem to reveal the interaction effect between the 
DW particles in the zigzag-like quantum spin chain. 

In addition, we would like to comment on the recently 
synthesized zigzag compound (N 2 H 5 )CuCl3, which is 
considered as a isotropic zigzag chain of the coupling con- 
stant J1/J2 — 0.25.— In this parameter region, the mag- 
netization curve does not have the 1/3 plateau. However, 
the compound belongs to the region of the even-odd os- 
cillation. Of course the bulk magnetization curve does 



not show such oscillating behavior, but we can expect 
the dual structure of the DW and its bound state in 
a magnetic field. According to the present theory, the 
magnetic excitation accompanying the ±1 change of the 
magnetization should be described by the combination 
of S z = ±1/3 DW and S z = ±2/3 DW. Recalling that 
ESR or NMR measurements is theoretically described by 
(S + S~) type correlation function, we can expect that the 
spectra of these experiments may resolve the DW exci- 
tation and its bound state. 

Finally we wish to stress that our approach is quite 
general and is applicable to the plateau states represented 
by the Ising bases accompanying the spontaneously sym- 
metry breaking of the translational invariance. The DW 
particles having the fractional S z is essentially originat- 
ing from the discreteness of the Ising bases and the de- 
generating plateau state; If there is a P-fold degenerat- 
ing plateau state, there is a possibility of a fractional 
DW excitation of S z = l/P- It is actually known that 
the S = 1 Ising-like zigzag XXZ chain exhibits various 
long-period state in a magnetic fields Our approach is 
useful for analysis of such a problem. In addition, re- 
cently, a J\f=2 super-symmetric ID lattice model is pro- 
posed, where the very similar ordered state and the ex- 
citation having the fractional U(l)-charge are explicitly 
constructed i^I It may be an interesting problem whether 
such mathematical model can be associated with the our 
explicit construction of the fractional S z excitation. 
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